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ABSTRACT 


fiteweiot-lolstoy (B-T) exact impulse solution of diffrac- 
tion by an infinite half-plane is compared to the usual 
HBeimholtz-Kirchhoff (H-K) integral formulation and to the exact 
mememmnuous wave (CW) solution of Macdonald. For backscatter 
meme -!| and H-K solutions are found to differ significantly, 
especially near the surface of the half-plane, where the B-T 
solution gives close agreement with experiment. For forward 
memecer the two exact solutions and experimental data are in 
agreement. B-T is found to agree well with measurements of 
diffraction by a barrier perpendicular to a rigid base. By 
memetaering source and source image in the base separately the 
Benmeept of “image of the source in the barrier" is found to be 
unnecessary. Use of the time domain form of the B-T solution 
in calculating the forward diffraction near a corner and be- 
fomdea thin strip is shown to give results which agree well 
with measured data. SOCOldamiwndI tanec tiOn errects are ob- 
served in the measurements of diffraction by a thin strip, a 


non-vertical barrier and a thick edge. 
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I. INTRODUCTION 


Many practical problems in acoustics are complicated by 
the presence of boundaries which influence the sound field 
in some way. When the direct path between source and receiver 
is not obstructed, these boundaries represent a perturbation 
@emene «cree field solution. If the boundaries are infinite 
in extent, the image or normal mode methods are often applied. 
Finite boundaries introduce the additional phenomenon of 
Mmeereractilon Or scattering, which is often treated as an 
madre tonal perturbation to the free-field plus reflection 
Solution. When the finite boundary is such that it blocks 
the direct path the problem is somewhat different. 

In this case, the only sound reaching the receiver is 
that which has been diffracted at the edge of the finite 
boundary (assuming the boundary to be rigid). In terms of 
mame rme noise control a finite boundary that blocks the 
Samect Path is often called a barrier or screen. When a 
Similar technique is applied to the reduction of unwanted 
sound in a shipboard sonar system, the boundary is generally 
@epled a baffle. The location and size of baffles on sub- 
marines and surface ships are severely constrained by the 
size of the ships themselves. As a result, the baffle is 
Sreeen close to the radiating source and the receivers close 


to the baffle. Care must be used in selecting a theoretical 





approach to solving the diffraction problem under these cir- 
cumstances since some of the techniques in general use involve 
Simplifying assumptions which mav not be valid under these 
meeorous conditions. 

Much work on the solution to the barrier problem may be 
monde in the noisé control literature. In this paper, two 
Soeeerioens uSe€d in noise control are examined, along with the 
little known solution by Biot and Tolstoy, and the theoreti- 
Gal results are compared. In addition, experimental results 
meereing to the application of theoretical solutions to 


finite barriers are presented. 








Mime LON CSAP PROACH 


A. GENERAL 

itesusual approach to solving the problem of diffraction 
Peeeoeplate Or barrier 1S tO approximate the barrier mathe- 
matically as an infinitely thin, opaque, half-plane. This 
problem was first solved by Sommerfeld [Ref. 1] for plane 
Somermmuous waves and then by Carslaw [Ref. 2] for spherical 
continuous waves. Macdonald's [Ref. 3] work for spherical 
continuous waves is often referenced as the rigorous solu- 
tion to diffraction by a wedge (and hence the half-plane). 
menos, Biot and Tolstoy [Ref. 4] published their exact 
impulse solution using the normal coordinate formulation. 
Despite the fact that the Biot-Tolstoy impulse solution is 
in closed form and contains only simple functions, many 
authors in the noise reduction field still consider Macdonald 
as the only exact solution [Ref. 5]. 

A second general approach to the problem is kirchhoff's 
approximate solution to the integral equation formulation 
of Huygens' Principle. This approach has been used exten- 
feemely in both optics and acoustics and can be shown to give 
reasonable agreement with observed results, within certain 
Mees , in spite of the fact that the basic simplifying 
assumptions are generally acknowledged to be incorrect 


fer, 6]. 





There are additional approaches to this problem that will 
femeee Explored here. A good summary of the various solutions 
memeene Context of the noise control barrier can be found in 
Reference 5. 

Rigorous derivation of the three theories will not be 
attempted, however the development of each will be summarized 


mimordcer that basic differences can be understood. 


Beeb lOl-TOLSTOY FORMULATION 
M.A. Biot and I. Tolstoy published their landmark paper 
Peewee +| On the normal coordinate approach to wave propagation 
in infinite media in 1957. Within this paper was an applica- 
tion of the normal coordinate approach in solving the problem 
of diffraction from an infinite wedge which forms the basis 
Mmemeeme theoretical results in this paper. The important 
features of their work in the present application are as 
follows: 
feeeeline SOlUtCLOn 1S exact and in closed form. 
b. The solution can be expressed in elementary functions. 
c. The use of an impulse source provides insight into 
the contributions from various edges and a means of 
Homabating the dittracted from reflected and direct 
components. 
The extension of their results to other types of sources is 
straightforward. 


immvesmormal coordinate approach is related to the normal 


mode description in acoustics which is most often encountered 


10 





in the solution of problems where the medium is bounded such 
as room acoustics, shallow water sound propagation, and 
Seeker enclosures. In these cases the normal coordinates 
are the allowable modes of vibration (normal modes) of the 
feecem. the overall response of the system is the super- 
Meereron of these discrete normal modes. Stated another way, 
mmmpeceeimal SOuUrce excites all modes at various amplitudes 
and phases and the solution then is the coherent sum of all 
meso individual excitations. The essence of the normal 
menaimate approach to the solution of general diffraction 
Meeeeems 1S the extension of the principle of superposition 
of normal modes to infinite media. As the boundaries of the 
Menclosure'’ are removed toward infinity the spacing between 
normal modes becomes progressively less, with the limit 
Beeeea COntinuous response spectrum with the boundaries at 
Memity., iIlf we now introduce a source which "turns on" at 
Mamet = 0 (a transient source), the solution contains only 
Pemwes traveling outward and no radiation condition need be 
imposed. The solution summarized here is taken primarily 
irom keLerence 7. 

In order to formulate the problem of diffraction by a 
wedge, a cylindrical coordinate system is chosen with the 
Origin at the wedge crest and the z-axis along the crest as 
shown in Figure 1. In this system the acoustic wave equation, 


in terms of displacement potential 9, 1S written as 


Isl 





+ Receiver ; 


Preure l. General Infinite :edge Diffraction Coordinate 


oystem. 














a*o , 1 3¢ 1 3*o , 3%9 1 3*o 
or? * y ar * r2 382 * 322” GF pez” 9 a 
Pyare paration of variables the harmonic solutions are 
vie eee 


Vv 


where the separation constants are related by 


Applying the boundary conditions for a rigid wedge, 
a _ , 
a a 


meemericting the amplitude to be finite at r = 0, and choosing 


feeeocation fOr Our point source (r = To? @ = BO Z4=- 0) eeinac 


femieeresult im Z-axis symmetry, the solution is written as 


i - 4,, (k r)cosv 9 cosy ze." (ese) 
n 
Vv, = a Cire sn 


W 
The solution can be separated into space-dependent and 


time dependent parts, 


on = Gav a (II -4a) 


<= 
il 


J, (<r) cosy Scosyz (II -4b) 
n 


The q,, are the normal coordinates which satisfy the 


differential equation 


: Qn 
2 = 7 
I, 7 8 dy UL, a) 





mere Qn is the generalized source function and i 1S pro- 
pemtronal to the kinetic energy in each mode. The series 
mepresented by (Iff-4) is known to be orthogonal and may be 


normalized by the condition 


+ co 


Ow 
i = 0 | alls, 


2 
dz | CW) aed as 


0 0 
Sieemedium 15 infinite in the r and z directions resulting 
in divergent integrals. These are handled by Tolstoy using 


a symbolic limiting procedure (discussed in detail in Re- 


Meemee, /, Section 8-2) resulting in the following expression 





w* TT 
ae 


w c* <dkdy oie) 


n 
Tt is at this point that we must choose an appropriate form 
fOr the generalized source, Qn 

The usual approach (and the one which we will ultimately 
use) is to choose the delta function or unit impulse source. 
Tt has two distinct advantages over other forms. First, it 
will result in a solution identical to that achieved via the 
Green's function approach, at least for those modes which are 
Maysitcally realizable. Second, it provides a simple building 
block with which the solution for other source forms can be 
constructed through the application of Duhamel's theorem. 
itemresulft 1s analogous to analysis of electronic filters: 


the response of the system to a unit impulse input is the 
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system impulse response. The system transfer function is 
the Fourier transform of this impulse response. 
Biot and Tolstoy choose the source such that the displace- 


Temempotential at a distance R is of the form 


a Ll 
2 4nR 


R 
cee Ue) einen) 
where l(t - Gy is the unit step function having the value 
zero for (t-*/c)<0 and one for i= ee 0) This represents 
Gie instantaneous injection of a unit volume at the source 
coordinates and is commonly used as a mathematical approxima- 


tion tO an explosion. Since the pressure and displacement 


are related by 





oO PAt2 Cis) 


the pressure at a distance R due to the source function of 
Equation (II -8) will be 


mee 
P “ GrR dt 


[6(t-*/C)] 

However, Medwin [Ref. 8] derives the solution for an impulse 
Pemnce directly from Biot-Tolstoy by simply taking the pres- 
sure to be proportional to the Biot-Tolstoy expression for 
Pemeicle velocity, since the two source functions differ by 
meeemrst derivative. Medwin shows that this is equivalent to 
assuming a point source of strength S (volume per time) which 


flows uniformly beginning at time t = 0. The acoustic 


pressure at range R due to this source is given by 


i 





fee ALC -t)Iv yd), (II -10) 


fememe A iS the bulk modulus of the fluid, and va 1s the space 

dependent solution for the free modes (Equation II-4). The 

mapaeian Of JW , 1n brackets, is evaluated at the location of 
n 


Mieesource. Substituting (]I1-4) into (II-10) and utilizing 


Meeeo, the differential Equation (II-5) canbe written as 


m2 
ie wg = = Cat) Jy, (<r Jcosv 6 Kdkdy Ce) 
W 


nh 





The solution to (II-11) represents the normal coordinates for 


the problen, 


Eos ute ZC * 
= ——— —— J (Kr_)cosv 8 Kdkdy t>0 Cliga 1.2) 
n me Tiere i O mo oe 
Sm@eeeituting into (Il-4a) and taking the time derivative, the 


Semeierdon can be written as 


me  (2C2 . a0 - 00 
a Te cosy 8 cosv 6. Be ey To! cosyz xX 
6lhC«OO 
Sin Wt cdedy aes) 


uy* 

Equation (11-13) represents an exact solution to the 
perturbation of a propagating wave by the presence of a 
rigid wedge, assuming Equation (II-10) as the form of the 
eommec. its present form is not particularly useful, how- 
ever the application of two known integral transforms to 
peertaon (11-13) reveals its simplicity and its close re- 


lationship to the physical problem in an intuitively 
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mor ying way. the application of these transforms is 
discussed next. 

Equation (II-13) contains integrals with respect to both 
the z-axis wavenumber y and the radial wavenumber kK. Con- 


Sidering the integration over y first, 


co co 
Sim We dy E 


wy 2 


i 
Soe yeahs) 9 Gy (TI -14) 


SOsSyY Z i 
(y 7+e 7)? 


cosyz 


A fh 


fleece relationship of the separation constants 
w? = cP (y 2+?) 
The integral on the RHS of Equation (II-14) represents a 


knownintegral transform, as follows [Ref. 7]: 


co 


Le 
: 2 orca L 
5 cosyz 212 ee ) dy = >a CmaboGeto-7-) | er>z 
caer lo ae 
) 
Sa OPe <7 iat 


imemexistance of two distinct forms of the solution, one of 
Which is zero, 1S a direct consequence of the choice of a 
transient source, i.e. if the source is turned on at t = 0 
there can be no effect at the receiver until a time equal to 
the source/receiver separation distance divided by the pro- 
pagation speed. In this case, z iS one component of that 
Separation vector. 

Using the result of (II-15) the remaining integral of 


(II -13) mav be written as 
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co 


s Pe 2 
i ey tr) spk cc mae) 


% 


]K ck (II -16) 
0 

Tas integral has a known solution which takes on three 

different forms depending upon the relationship of the 

variables r, r., and (c-t-= 22) 2, Tolstoy [Ref. 7] defines 

Mmrese three regions in terms of the propagation time of a 

pulse emitted by the source at t = 0. The time of arrival 


(at the receiver) of the earlist possible direct pulse is 


2 


—_ 


eeeven by L/ 


Be ((r-r,)? + 27] 


1 


Ale 


The earliest time of arrival of the pulse which has traveled 
from source to receiver via the crest of the wedge is 


pee 2 2 11/2 
cy aCe Ty + 2 | 


With these new variables, the three regions of the solution 


are given by 


ir , p<: I = 0 
. qo: 
-1 
ll. ect i - Gs sm Sink) cosv,% 
r? + Es peat 0<X <n 
MM = arcos ee 
° is) 
for | T <t i (aes ila Y)" sinty mye Yn 
0 nN 0 n 
c*t? _ (eras - ae i ne) 
Wee Bane cos i ITF 


CLL aes) 
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micise tEnree regions have a strong intuitive foundation. 
ieease |. it is obvious that a pulse transmitted by the 
weece Cannot be received until a finite time has elapsed, 
that time being equal to the separation divided by the pro- 
Pagation speed. Case II. represents that time during which 
reflections from the face of the wedge are received. Tolstoy 
[Ref. 7] shows that when this expression for rT. LomcUbSt hENted 
into Equation (II-13), the solution takes the form of images 
Seemerie SOUrCe in the reflecting wall. Finally, Case III. is 
the diffracted wave. It is characterized by an arrival time 
Weeene is jater than the last possible image source arrival, 
as seen from Figure l. 

inewsame results can also be interpreted geometrically, 
as discussed by Watson [Ref. 9]. If we rename the variables 
r, r,> and (c2t2-22)2 iti oymas ob. CC, amd a; respectively, 
MeesemebOInts out that the three regions of the solution can 
Mmemmaescribed in terms of the geometry of Figure 2 (z 1s taken 
memme zero £0r Simplicity). The solution of Case II. applies 
whenever the value of side a is such that a triangle can be 
memmed having sides a, b, c. It is obvious that this will 
be the case whenever a has the values 

foac )<a<(btc) 


and, from the law of cosines, the angle A will be given by 
eae AZ 
ZDE 


where A can take on values between o and 7. The comparison 


A = arccos (Ti -20) 


Pmeeecauatzion (11-18) is immediate. In Case III., a 1s 


Le 





A 
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—-> 
O 
U 
ae 
Sc 
e 
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i 
= 
ic 
ww 
at 
° 
= 
ce 
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an 
da 
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{ Eg 
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LL 
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e 
O 
m 
e 
ie 
ad 
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is 
zl 
a 
LY 
S 
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Meter than (btc), the angle A becomes imaginary, and the 
ewepeession corresponding to (11-20) must come from hyperbolic 
trigonometry: 


Bee oe 


A = argcosh te 


ube vie 


Memean be seen from this discussion that the Biot-Tolstoy 
formulation in the time domain decomposes into direct/reflec- 
ted and diffracted components in an intuitively satisfying 
way. This iS apparently due to the use of the normal coor- 
aoeaee) tOrmulation, which considers only modes which are 
Mieestecally realizable, and the choice of the transient form 
Pemmemen source. However, because of this distinct division 
between direct and diffracted solutions, considerable care 
fiiewbe exercised when working with a geometry where the 
arrival times of the direct and diffracted waves are close, 
such as near the boundary between the "illuminated" region 
and the geometric shadow. 

imomarrive at the final solution for the diffracted wave, 
m@emresult of Equation (11-19) is applied to the general 
solution (II -15) as follows: 
ve \ 


a> . “(rr Sent cosv.§ cosy 6 sinv me n 
O iy 0 rh n 


3 ie 
W 


— 
it 1 8 


ae 0 Q 
y) == Chic 72\) 
n Ow 

Expressing the sum in terms of exponentials, regrouping, and 

Bemelecting in conjugate pairs, the solution can be written 


meexact, closed form as 


ok 





ESS Tl Oe 


ae 1 1 +O+ 
wy sin (4 (74926 )) 


QQ 
©- 
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(rr. Sinh Y) 


Q? 
ct 


2 Pai 


Sw 
1-2e "Weos(7-(m#048)) +e . 
Ww 


eS 
479 

W 

(11-23a) 


c*t? ” (ies an a i no 


Y = argcosh (I1-23b) 


oi 
O 


memes the term (7™+8+9,) indicates that the entire expression 
meen cne Curly brackets is the sum of four terms, each of 
poem corresponds to a unique combination of (7+6+6,). 
Recalling the discussion concerning the forms of the 
impulse source, Equation (I1I-23) can now be used in place of 
Mmipmeeecond partial derivative in (11-8) and the solution of 


the diffracted wave in terms of pressure can be written as 
Pe 
p(t) = - eee 5 sinh Oy e ow [2 } 


g 
Ww cee 





where Y is given by Equation (I1-23b) and 8 represents the 
Berm in curly brackets in Equation (II -23a). 

meeeron (11-24) represents the diffracted pressure as a 
mmetion of time due to a unit impulse source. To within a 
constant representing the source strength, this is analogous 
to the impulse response, h(t}, of a general linear system. 
This impulse response is often Fourier transformed to obtain 
m@emtranster function H{w) and examined in the frequency 


domain. Since there is no known analytical transform of 
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Seeeoguation ([1-24), the numerical technique of computing the 
ieese response in discrete time steps and then transforming 
Demmethe Fast Fourier Transform (FFT) is employed. To avoid 
the problem of infinite values of h(t) for the earliest 
diffracted arrival, Medwin uses an approach which computes 
@iemitirst time point by a numerical integration. All of the 
Ieee )Oistoy theoretical results presented in this paper are 


@eremlated in this way. 


Gee te oMAOLTZ-KI RCHHOFF FORMULATION 

The Helmholtz-Kirchhoff formulation is a mathematical 
memeemenmt Of the heuristic description of wave propagation 
due to Huygens. Simply stated, a source at some distance 
from the surface (in this case a plate which has one straight 
edge and extends to infinity in three directions) insonifies 
all points on the surface S. Each small area element, ds, 
Smmene surtace then acts as apoint source of spherical waves. 
Peesome Observation point 0, the acoustic field is due to the 
incoming waves from the source directly plus the sum of the 
contributions from all of the point sources on the surface. 
Neglecting the effect of the direct incoming waves at 9, the 
field can be determined by integrating over the surface. To 
@emstruct an exact solution in this manner one would have 
to know the amplitude and phase of each point "source" on 
Mi@ewstrface. This not possible, in general, and therefore 


Ememaarchhotft approximation is used. The brief development 





meee herein follows that in Clay and Medwin [Ref. 10]. A 
meeeled derivation can be found in most optics texts and 
in a book by Baker and Copson [Ref. 6]. 


From Green's theorem, assuming a source of time dependence 


a, Sime tiled at the observation point Q can be expressed 
as 
= eels etre 
1 me: e 3U 
U(Q) flue - eda 
4m an R R on S (II -25) 


5 


where R is the distance from each element ds to the observa- 
-ikR 


R 


That the integrand is evaluated on the surface is denoted 


meme point Q and Powtmemineewspace Green's function. 





Meeeere SUbSCript s on the brackets. This known as the 
Helmholtz-Kirchhoff integral. It describes the relationship 
Of the field at 2 to the field on S, but the problem of 
evaluating the field on S is not yet resolved. 

Equation (II -25) can be used to evaluate the diffracted 
wave field on either side of the surface. Our primary in- 
terest is the diffracted field when source and receiver are 
Somopposite sides of the surface, however, when they are on 
the same side (herein referred to as backscatter) the result- 
img expressions can be somewhat simplified. Since the purpose 
of this development is to compare the results with those of 
the normal coordinate approach, the simpler case will be 
used. For a complete treatment of this Helmholtz-Kirchhoff 


itegral equation approach, see Baker and Copson [Ref. 6]. 





et U. Pomcectmiod as the imcident wave field at the loca- 
tion of the surface, with the surface removed, then the 
meemeaty U within the integral can be approximated by 

U = MU . (II -26) 
femme “ 1S the plane wave reflection coefficient. From 


Somerton (I1-26) the normal derivative can be written as 


JU dU 
on : on (11-27) 


Peeerdudation (II -26), the restriction has been imposed that 
Mm@emonly wave field that can exist on the insonified surface 
mmerme Darrier is that due to reflection of the incident 
wave field (i.e., there can be no scattered field on the 
Surface). The consequence of this assumption will be seen 
mmemecme diffracted field near the surface is examined. 

The additional assumption that is required to make the 
Meeceral in (11-25) tractable concerns the reflection co- 
eeereient M. Kirchhoff's assumption considers it to be 
equal to the plane wave reflection coefficient for an infi- 
ieee plane interface. Here the surface is assumed to be 
Memrectiy rigid and M= 1. Equation (JI -25) can now be 


written as 
ac 
Bl 3 elm 
mea) = an | nis TS 





(II -28) 
When the integral in Equation (II-28) is taken over the 


miierace of the half-plane, the U(Q) represents the backscat - 


tered wave field, If, on the other hand, the integral in 


Pouation (11-28) is taken over the portion of the plane not 


Zs 





Seeupi1ed by the surface (generally referred 
mao ptics literature), U(Q) would represent 
meeecered torward into the shadow region on 


fewer -Hlane opposite the source. (Note that 


to as the aperture 
the wave field 
tiie aside "Or. tie 


this approach 


foc also require a change of sign due to Equation (II-27). 


D. MACDONALD'S FORMULATION 


The final theoretical development to be considered here 


memmeeiae Of Macdonald [Ref. 3]. His solution, like Biot- 


Memeeroy, 1S an exact solution for the infinite rigid wedge, 


except that Macdonald assumes a continuous wave (cw) source. 


Peaerivation of his solution will not be attempted, instead, 


the results of Kawai, Fujimoto andItow [Ref. 


11) will be 


Moeae the calculations in their paper are based on an ap- 


proximate form of Macdonald's solution given by Bowman and 


Pemmenr |Reft. 12] but only the form of the exact solution is 


Meesented below. 


From Reference 11 the exact solution can be written in 


terms of the velocity potential as 


< 
il 


vce.) + V(-8,) 


‘ ©0 (1) 
feat ls + kR) 


H, (T?* + kR) 


v = ik as ols | dT 
meee 1 oie. on; 
as Pera é meet 2kR)+ é 
a 


-T, pe 


where 
ie 


ees oen(t - 8 -6) (KR, - kR) 
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(11-29) 
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: 
Tae = sgn(t ~ @ — 6.) (kR, - kr*)}/? (II -30b) 


ey et 
erm) = 1. 1 X<0 
2 2 2 ie 
hoe (rr + ea 32 eC cos(é oa) (II -3la) 
R* = (r? + 2 i>. 2 Lf 
iG ie Za Ze cos(6 + 8 I) (II -31b) 
, . ey 2 
meme (r er je + 2°} Gils ie) 
k = w/c 
Tr, To» 2; oe 8 dome rimnedmain Ficune ~ I 


and ee!) ete rince Ondemmneanke!l fTumction of the first 
Keer cl. 

[ijemoreniticant feature of this solution, first pointed 
out by Macdonald [Ref. 3], is that it consists of the sum of 
two velocity potentials, V(9.) and V(-8,). These are iden- 
fjered in the noise control literature as the source and 
mem tmage Of the source in the barrier", respectively. This 
mpenmcant interpretation of the solution will be discussed 


fea later section. 





Cit wecoren nt oON OF SOLUTIONS 


Peet OL-TOLSTOY AND HELMHOLTZ -KIRCHHOFF 

Peort-1olstoy (B-T) and Helmholtz-Kirchhoff (H-K) solu- 
Enons Willi be compared by choosing a geometry in which both 
Can be expressed in closed form. The backscatter case is 
chosen wherein source and receiver are in the same location. 
For this case the time domain solution for the H-K formulation 
as derived from Equation (II-28), once again following the 
development in Clay and Medwin [{Ref. 10] with some minor 
changes in notation. 

Taking source and receiver to be in the same location 
(backscatter only) and assuming an impulse source, U_ in 


> 


Equation (II -28) can be written as 
e ikR 
U. = = GLE) 


Making this substitution and transforming the surface 
integral into an integral along the edge [Ref. 10, pg. 323] 


Equation II-Z8) can be written as follows: 


S 
aa 


U(Q) = JJ 4 jac 


Per 2) 





[emmeny, ©, and R are defined as shown in Figure 3. If the 
fees) line through 2 intersects the rigid half-plane, ene. 


integrand of equation ({III-2) is represented bv 


Me ms 
me hUTh hE] 
co iG 
(Case [) 
[If not, the integrand becomes 
r 
[ 
(Case II} 


(iwemeetie inverse Fourier transform of (111-2), and denoting 


the time domain solution by a small u, 





= 
i ee oe 
u(Q) = = | fs —__— aS df da 
Re 
-@O-T CO 
a 
Guinean) 
For Case I this becomes 
coe 2) Ge at) io (* ji (wt-2kr) 
4 OTSA a 
u(Q) [ ny ie ae ‘I ea er 
col 
2 
Cir ) 


Where 6 represents the Dirac delta function, a consequence 
of the choice of the Green's function of the Helmholtz equa- 
tion for the incident wave field. 

Mme rtfirst bracket on the right hand side of (III -4) 


represents the reflection from one nalf of an infinite plane. 
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pe 


“d 


Figure 4-6 Coordinate System for Backscatter from an 
Infinite kdge. 





u(Q) 


u(Q) = Fey D(t) 


Mmessecond bracket contains the specular reflection from the 
other half of the plane but reduced by the integral expres- 
Meomerepresenting the fact that the plate is finite. This 

meme the torm of an integral along the edge. Following Clay 
and Medwin [10], this is called the boundary wave and defined 


as follows: 


2 1(wt-2kr) 
Dt ) = ef / ———— dadf 
- © ieee 


ul 
2 


r| a 


(PiEss ) 


M@eoe if results in a similar expression but without the 
reflection. 
In terms of Equation (III-5) the following expressions 


mm@emmritten for the time domain impulse solution: 
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Bee 
age) 
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Taking advantage of the time domain form of the solution 
amd Separating the reflected and diffracted (boundary wave) 


components of the wave field results in 


(-) Case I 


nee y. 
ot ee, Pe (+) Case II 





joes can be written in terms of the geometry of Figure 3 as 


follows (from reference 10, section 10.2.5) 


ol 





ak 


y x 
mo) = : wees = SLT } 
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(III -8a) 
= 0 iE ho (III -8b) 
where 
2x, 
Be = meee eae COSmor == Tp COS oF 
~ CT» 
” 2 
2R 
4 = ——— 
E 
Mmiomarly, through trigonometric substitutions, 
2 hoesin 29 
a) 3 { - 2 2 ey 2 


Miemalttracted pressure at Q as a function of time is given 
by 

Pqg(t, Q) = - PsRiu(Q) meena 
where Ps 1s the source pressure spectral density and se 1S 
the reference distance. Taking the product PR. = | Pascal 


Meter, the diffracted pressure at Q, as a function of time, 


Simeone written as 


Pg(t, Q) ce’ u(Q) 
2 bor Sin. Ze 


= - _ ere sere es 


epee eE a e  nea 


Gi eteesinl 


iy 
tO 





miemcOnresponding expression for the Biot-Tolstoy for- 


mulation can be obtained from equation (II-22). 








For the 
thin plate ce Se 
ny 
ae. C (er. sinh Y) © cos ne EOS a Sa) ae ea 
a 272 =« n " : Cum 12) 


Wemking with the sum and noting that 6 = 8, (source and re- 


ceiver coincident), 


ae in§g -inée ans 
° 7 = e a 
5S = Ecos 282 sins e "~ = rs (gle paint Ci 2 ae e 7 


a a 2 Z 


pultaplying terms and observing that sin nv is zero for n 


Z 
even and +l for n odd, the sum can be evaluated as 
ees oY 
S = 5 ( e a 2 as e [ee a! ak 
+ z ee ee OT se 
+ + Mer - e > 2 + g Ue. Le ) 
ee i 1 1 ' i 
7 sech ano sech uu, + 3 sech uz; Re (5) >0 
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ay = ig Sec 
Uo = 5 x 16 
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After considerable manipulation the sum can be written as 


t E cosh? COcmUmENeGsmisy + COS 26 
s =f [——=—________ 
(cosh Y + cos 280) coshs (III -14) 


-_ 


From the geometry of Figure 3 the following can be defined 
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iiemvarious trigonometric quantities in equation (III-4) can 


menace expressed in terms of these as follows: 


ao (= - Lo 
to* 


oa 


mes 260 = |] - 2 sin*6 =1- 2 


ror rT = re eit ronme( fl -Zob) can De Written as 


fe= argcosh (V) 
mee co tO 7 ae 
maere V 7 - 1) = (2 = 


Using the half angle relations for hyperbolic functions 


- 1) 





leads to 
, 
cosh 5 - Scve1)+/%) - ae 
pumtlarly , 
Pm = (v2- 1)2/2 
~~ Se 2 ely 
7 Le Gee, ) 
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macetinally, to within a factor of 2, 
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betting PRG be unity, as before, results in an expression 
for the diffracted pressure wave in the time domain which 
Saaepe Compared directly to the H-K result. 


iot-lOolstoy: 


2 . Beats 
a ee — 
iGuameein teem T, ) Gee 6) 
memmmoltz-Kirchhoff: 
Z 
p(t) ag meee We oa a t>T, 
@ites nice 1/2 CU i ) 


memaeion (Jil-16) and (111-17) are similar in overall 
mOrm but differ in several important ways. When ce is much 


Peeater than to: the time dependence of both expression is 


X 


determined by the Geary) factor. (for t, approaching 
TOs the behavior is more complex and is treated in detail 
maeexppendix A). 

Since the impulse source has an infinite amplitude, the 
leading edge of the diffracted pressure signal (tec) must 
also be infinite. After this infinite leading edge, the 
time behavior can be described as follows: 

2 


At = t - To e- = hee 7 he he 


Memecne ‘early diffracted signal, ac and 
-1/2 
: Sole 

Ct? = 02)? = (2tt, - 2002)? = Cate (t-to)) 1! 2= (27, At) 
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Mimeminitial diffracted time signal decays as At * for both 
Seeeerons. Since the amplitude is highest in the early time 
pomeenon Of the signal, this At? dominates the Fourier trans- 
form of the signal and the frequency dependance of both 


ae 
mamicions takes the form f °* 


Wiwhen tS -Cnardeteristic of many 
diffraction problems. As time increases, however, the H-K 
solution approaches the form t" while the B-T form approaches 
a 

Both solutions have been implemented on the IBM computer 
meweomeat NPS. The diffracted time signals for a representa- 
tive geometry are shown in Figure 4. Care should be exer- 
=erodman comparing the magnitude of the first point (t = 0) 
Since the "infinite" values here were handled differently in 
the two solutions. The approach used in the B-T computation 
aSyene numerical integration technique discussed earlier 
while the apvroach for the H-K results are taken directly 
from Clay and Medwin [Ref. 10]. 

iiespehavior ofthe two solutions as a function of 6 is 
also fundamentally different. Examination of the numerators 
of Equations (III -16) and (III-17) over the region of 9 from 
0 to 180° (the problem of the thin plate being symmetrical 
with respect to 180°), reveals three apparent zero values 
momeecme H-K solution (8, 90°, 180°} vice two for B-T (90°, 
meee Actually, the zero in the numerator at 90° is offset, 


in the limit, by the term (t*-t *) in the denominator which 
0 


Mees tO infinity at t = to: Physically, this means that 
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Mem SOurce and receiver are directly over the edge, perpen- 
ieeemrar to the plate, the solution consists of reflection 
only. In order to show the detailed @ dependence, the time 
domain solution has been computed for a variety of geometries 
and each time solution transformed to the frequency domain. 
These results are plotted versus 9 for several frequencies 
in Figures 5 and 6. The amplitudes are in dB with the 
Memerence being twice the range from source to edge. The 
differences between the two solutions are significant, with 
the only region of agreement being near 90° where both solu- 
Mm@omonappear to approach the value of -6dB. This is the 
expected result and represents one half of the pressure 
maigemeude that would result from reflection by an infinite 
Mmeaews ACtually, the H-K formulation produces the correct 
result at 8 = 90° while the B-T solution appears to fail in 
the immediate vicinity of 90°. As 8 approaches 90°, the 
travel time of the reflection ty approaches the travel time 
of the leading edge of the diffracted pulse oe 

In this region, equivalent to approaching the shadow 
boundary in the forward scattering case, the B-T solution is 
approaching the delta function form and the approximation 
Meead to calculate the first time point may not be valid 
(see Appendix A). The B-T diffraction solution can only be 
Metea for © approaching 90° since at 9 = 90° the separate 


meeiection solution would apply (Equations II-18 and II-19). 
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Peete Source and receiver are moved away from 90° the 
feereoreement between the solutions is dramatic. It is in- 
Mueeeevely correct that both solutions go to zero as @ 
meeroaches 180°, since the plate is infinitely thin and 
would not be "sensed'' by a wave from this direction. The 
way in which the solutions approach this zero value is 
mueece different. 

Miewmost significant disagreement between the solutions 
mempsecl when source and receiver are close to the plate (98 
approaching 0). The H-K solution approaches zero near the 
Metees (it is svmmetrical, except for sign, about 90° due to 
eiemomm 26 term) while B-T approaches a constant value for 
Peemeiregquency. This zero in the H-K solution appears as a 
eemeegquence of the simplifying assumption stated in Equation 
Mieco), where the total pressure field in the plane of the 
meteor o assumed to be the value of the incoming wave field 
evaluated at that location (M = 1 assumed). This "boundary 
Gonmd1tion'' therefore forces the diffracted pressure field 
mummeemzero anywhere in the plane of the plate. It will be 
Shown in the experimental portion of this paper that the 
Pec Oolhution is correct in this region near the plate. 
Peeeomen not shown in Figures 5 and 6, both solutions exhibit 
m@empnoper Sign behavior, i.e., the direct and diffracted 
flises are Opposite in sign in the region 0 to 90° but have 


the same sign for 90 to 180°. 


4) 








ieescummary, the B-T and H-K solutions simplified for the 
Meer Ot Coincident source and receiver, are dramatically 
different. They approach agreement, in general, only in the 
meen near, but not at, 6 = 90°. 

The final consideration before leaving the topic of 
fiescatter iS an examination of the form of the B-T solution. 
Multiplying out the numerater and regrouping terms, Equation 


Wemeet6) can be written as 


Zan cos 6 


1 

ect) mea | 5.0 iar (eee ar Sey eee 

aa?! Pome) (C7 ty) Bee ot) 
erie) 


When written in this way, the solution is seen to have an 
Mieemeesting Characteristic. The solution is the sum of two 
terms, one which takes into account the overall geometry of 
source/receiver and plate (the first term in III-18) and a 
term which depends on the distance from the edge of the plate 
to the source/receiver. This seems to be a characteristic 

meee xact SOlutions of the problem of diffraction by a half- 
plane and is directly related to the "image in the barrier'' 
first mentioned by Macdonald [Ref. 3]. This "image" approach 


homer be discussed in some detail in the next section. 


Beeeeol OT-TOLSTOY AND MACDONALD 

In comparing the B-T and Macdonald half-plane diffraction 
Seutions, an approach similar to that of the previous sec- 
Mion will be used. However, the solutions will not be limited 


to backscatter but will be kept general and the emphasis will 
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Memon = LtOrward scatter", i.e., diffraction of sound into the 
geometric shadow region (source and receiver on opposite sides 
Sees half-plane). The purpose of the comparison is twofold: 
to show that computations based on the two "exact" solutions, 
one impulse and one CW do, in fact, produce the same results, 
and to gain some insight into the use of the "image in the 
momeerer'. AS discussed in section II, the Macdonald's 
solution form and computations are from Kawai's excellent 
paper [Ref. 11). 

The comparison between B-T and Macdonald computational 
memes Frequires some manipulation of geometric quantities. 
ijewat s results are converted to the notation of this paper 
for consistency, however, this results in only 2-3 points 
Der ftrequency spectrum. To alleviate this problem somewhat, 
tae geometry is selected in such a way that several different 
B-T runs result in basically the same values. Figure 7 
Shows this geometry and the comparison between B-T and 
Macdonald results. The reference pressure is the free-field 
pressure that would be measured at a distance of r + er from 
the source with the plate removed. Data referenced in this 
manner are labeled "dB re free field at 2r". This combination 
of geometry and reference results in B-T theoretical values 
mame vary only one dB over the range 0 = 330° to 9 = 360°. 
Kawai's computations are then converted from his normalized 
Meesentation {[Ref. 11, Figure 3] to the format of Figure 6 
° increments of @ between 330° and 360° inclusive and 


in 5 
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plotted in Figure 7 as discrete points. The agreement between 
the two solutions is seen to be very good. The scatter in the 
KaWai points is almost certainly attributable to the process 
Seeoeading these values from the graphs of Ref. 11. Additional 
comparisons were made over a wide range of different geometries 
Beeaeethne B-7T theoretical results of Bremhorst [Ref. 12] with 
Seeeement Similar to that shown in Figure 7 observed in all 
cases. 

Mmeesnould be noted, however, that the comparisons were 
not made for the general case of high frequency, large range, 
near the shadow boundary (Ry large and N small in the presen- 
tation of Kawai). 

Kawai [Ref. 11] actually presents two sets of computations 
of Macdonald's theory; one for the "real" source only (refer- 


mea tO as ATT, Jegand one for the "real" plus “‘image' source 
O 


(ATT). This terminology appears to be traceable directly to 
Macdonald [Ref. 3]. In his solution he separates a complex 
Mreeotal {Ref. 3, p. 422] into two integrals, performs the 
integrations separately, and recombines the two into an over- 
mesoOlution for velocity potential [Ref. 3, p. 425]. Bowman 
megmoenior's approximate form of this overall solution [Ref. 
Mmeis given as equations (11-29), (11-30), and (II-31). 
Mmiemcomparison of B-T and Macdonald in Figure 7, is based on 
mmecaonald's complete solution, ATT. Indeed, if the complete 
meecion of the half-plane problem is the sum of these two 


terms, it is not Clear why one would want to consider them 
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Meeerately. Yet this is precisely what is done by many authors 
Meeene noise reduction field [cf. Ref. 5]. 

In order to examine the "real" and "image" terms in Mac- 
womata' S SOlution more closely, equations (II-30) and (II-31) 


meemecpeated here with a slight change in the way 6 and 9 


O 
are expressed; 
Tp = sgn(m - (@ — 8 J) (kRi - ka) l/? (II -30a) 
= alee 
ee Ooh tm — (9 + 8. J) (KR - kR*) (II -30b) 
eee te ry tz? Cr were cos(e—9))1/? (II-5la) 
[ee (r°> + ro* + 22 - err. cos(9+8.))1/? (11 -51b) 
The primed values represent the ''image'' source terms. The 


real source terms are a function of the difference between 
8 and a while the "image'' terms are a function of their sum. 
iitemsteniticance of this is shown graphically in Figure 8. 


Figures 8(a) and 8(b) show two different geometries, both of 


which have the same value of the quantity 9 — 8 Simec 
(@ — 9 A) is only form in which the quantities @ and o enter 


Mmemeesource portion of the solution, it is easily seen that 
the real source solution is independent of the orientation 
of the plate. It depends only on the locations of source 
and receiver relative to the edge. However, the geometries 
of Figures 8({a) and 8(b) are different in terms of (9 + 9) - 
Tt is this term which retains the information concerning the 


Bytentation of the plate (half-plane) relative to source and 
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Pigure 8. Geometry of the Image of the Source in the 
paeimbe i. 
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Meeedver. Since changing the sign of 8, 1s the same, geome- 
Mercally, aS moving the source to the opposite side of the 
barrier, one can easily see why the second term of the solu- 
@rome 1S associated with the “image of the source in the 
Mmeeraer’. This interpretation is shown in Figures 8(c) and 
ered ) . 

Mio interpretation 1S not unique to Macdonald's solution. 
Recalling equation (III-18), is was shown that the B-T 
pemtion, for the case of coincident source and receiver, 
could be written as the sum of two terms, one which depended 
only on the orientation of the edge and one which also con- 
Sidered the location of the plate. Indeed, the general form 
Mueeee B-1 solution given by equation (II-24) could be 
written as the sum of a term that is a function of 6 — 8 
Mmemanother that 1s a function of 6 + 8 This would be done 


in the following way: 
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where each 8 term is actually the sum of two terms, one with 
meemcerm in brackets evaluated using the (+) sign and the 
@emer using the (-) sign. There appears to be no advantage 
Memwiating the B-T solution in this way except to show the 
@eemetrical Similarity to Macdonald's result. Indeed, Biot 
Gee OlStoy specifically identify both terms as diffraction 
Memerer lection. True reflections have separately been iden- 
meted bY B-T as occurring before the diffraction; an inter- 
Meeeerion that is direct and obvious in their impulse solu- 
tion, but which is obscured in Macdonald's CW solution. 

In summary, the B-T and Macdonald exact solutions have 
meen £ound to be in very close agreement for all of the cases 
examined. Real and image sources have been shown to be a 
merit not particularly useful, interpretation of both 


Somutions. 





iv> EXPERIMENTAL APPROACH 


fee GENERAL 

mrenmeasurements reported herein were made in the NPS 
anechoic chamber utilizing pulse techniques. The anechoic 
chamber provided the low noise, low reverberation environment 
necessary to make precise amplitude measurements at low source 
Memes. ithe use of a pulsed source provided a signal that 
was both wideband and deterministic and made it possible to 
Mammen ol the effects of extraneous scattering by judicious 
mmomee Of the sampling window. Detailed descriptions of the 
anechoic chamber and the data acquisition hardware may be 
found in Reference 12, and will not be repeated here, except 


Pmememecessary to explain the results. 


Pee ALA ACQUISITION 

A block diagram of the overall experimental data acquisi- 
f[eemeprocess is shown in Figure 9. Overall timing of the 
acquisition process is controlled by the programmable timing 
simulator. A typical sequence is as follows. The waveform 
generator is triggered by the timing simulator, sending a 
single replica of the preprogrammed pulse to the power ampli- 
meer, This pulse is amplified to approximately 50 volts 
peak, added to the 150 volt D.C. bias voltage, and applied 


to the source (source and receiver characteristics will be 
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Mmscussed separately). The transmitted signal is received by 
the % inch condenser microphone (BtK 4133) and sent to the 
Signal conditioning equipment. Here the signal is amplified, 
mepcered to remove low frequency noise and to prevent alias- 
ing in the Sampling process, further amplified and digitized 
(12 bits). Here the timing simulator determines the start 
time and duration of the sampling of each pulse. By selecting 
the sampling frequency and duration, the number of samples 
available to the FFT can be determined. In practice, the 
number of samples is set slightly higher than the number to 
De transformed since the computer begins with the first 
Sample and disregards those coming after the specified number 
have been acquired. This process is repeated for each pulse 
With a pause of approximately 80 milliseconds between pulses. 
The computer accumulates the samples by averaging up to 
10,000 samples in the time domain and then computes the 
Fourier transform of this averaged waveform using a standard 
Sertware FFT algorithm. Since the FFT must be taken only 
once, the system runs in real-time, which allows a very 
large number of pulses to be averaged in a reasonable amount 
St time. Since this averaging is coherent, the signal-to- 
memese ratio is proportional to N, the number of pulses, 
father than the usual N2, The averaged time waveform or 


mereresults are stored on a floppy disk for later use. 
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meee SOUND SOURCE 

The ideal sound source for these experiments would be a 
Meant SOurce capable of high acoustic output over very short 
time durations, with zero output at all other times. These 
characteristics can only be approached with a real source. 
The obvious choice of a spark source is rejected in order to 
avoid finite amplitude effects. Fortunately, the extremely 
low noise levels of the anechoic chamber plus coherent 
averaging makes accurate measurements possible with a source 
of relatively low intensity. The search for a useful source 
of small dimensions and narrow time response has been a con- 
tinuing one at NPS[cf. Ref. 14]. The source chosen for the 
experimental work in this paper were the B & K ¥4 inch and l 
milemeconaenser microphones. The use of these microphones as 
sound sources is described in Reference 12 and in a B § K 
technical pamphlet [Ref. 15]. In general, the % inch source 
covered the higher frequencies (6-40. kHz) while the 1 inch 
source was used to provide additional low frequency coverage 
Beem tO approximately 1.5 kHz). The 1 inch source is 
meetiatiy capable of providing sufficient output over the 
entire frequency range but directivity does become signifi- 
@ant at the higher frequencies. 

Typical source waveform and frequency spectra are shown 
meerieures 10 through 13. In Figure 10, the waveform at 
the output of the power amplifier is the top curve with the 


received waveform shown below. The received waveform (and 
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meoclated frequency spectrum) is actually due to a combina- 
Mmeoneot sSOurce and receiver response, however, the relative 
Contributions of each are not particularly important since 
meee odltsS are normalized by an appropriate reference mea- 
Surement using the same source and receiver. The time dura- 
mueneecor ringing") of the pulse is important since all of 
miem@esired diffracted pulse must be sampled before the 
unwanted arrivals from the other edges of plate. The longer 
EaewaGOuStic pulse, the larger the specimen must be to 
Bewemeare the pulses. This ringing is especially important 
Mieasuring backscatter since the low amplitude diffracted 
pulse will arrive after the direct/reflected pulse which is 
much higher in amplitude. Figure 11 shows the energy 
“eeemral density of acoustic pulse in Figure 10. 

mreures 12 and 13 present the same data as Figures 10 
enaetl, with the same receiver, but using the 1 inch source. 
In this case, the received pulse was low pass filtered at 
nieeieeanm order to enhance the low frequency portion of the 
SOc mr UM. 

Mice an acceptable source and receiver were selected, a 
mounting system had to be devised which allowed the micophones 
to be positioned accurately yet present a minimum of 
@dditional scattering surface. The approach used is shown 
Mierieure 14. The diffracting object (in this case a thick 
plate) was suspended with nylon fishing line from a frame 


attached to the ceiling of the achechoic chamber. Heavier 
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Figure 14. Typical Anechoic Chamber Measurement Apvaratus. 





@emectS were also supported from below. The source and re- 
Meaver were each attached to a pair of thin stainless steel 
fees mounted firmly to the ceiling frame and resiliently 

to the wire mesh floor of the chamber. Source and receiver 
moemeece positioned in the horizontal by moving the ceiling and 
floor attachments and vertically by moving up or down the 
wires. Additional support of the signal cables was provided, 


memetvecesSary, with nylon line from the ceiling frame. 


ieee LATE CONSIDERATIONS 

The plate used in the experiment should be as thin and 
vet, as rigid and non-transparent as possible. In addition, 
1t was desirable to keep the overall weight of the plate low 
@emeacilitate handling and suspension within the anechoic 
eaemioer., Aluminum plate .475 cm thick was chosen as a 
reasonable compromise. The transmission loss, calculated 
using the approach of Reference 12, was greater than 40 dB 
at 1.5 kHz, increasing with frequency. This was considered 
adequate since no thin plate measurements were planned 
where the direct path through the plate was to be substan- 
tially less than the diffracted path. For the thick plate 
tests the direct path would be short but two plates would be 
Meed and the resulting attenuation would certainly be adequate. 
These estimates were verified by examining the received sig- 
mmlmorior to the arrival of thefirst diffraction to ensure 


there was no significant direct path arrival. 
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Deere o iE RTE PERIMENTAL RESULTS 


In Figures 5 and 6 it was shown that the B-T and H-K 
memectons give dramatically different results for backscatter 
mear the surface of the rigid half-plane. A simple experi- 
ment was performed to determine which solution is correct. 
The one inch B+K microphone was used as the source with the 
one-half inch B+K microphone used as the receiver. They were 
taped together and positioned at a distance of 25 cm from the 
edge at various angles as shown in Figure 15. Although the 
source and receiver acoustic centers are separated bv apprcxi- 
ieeeely 2.5 cm they will be assumed to be coincident. 

ire transmitted waveform was the halfJtriangle discussed 
Carlier but the received signal was low pass filtered at 
60. kHz vice the 14. kHz shown in Figures 12 and 15, to 
achieve a higher frequency response. Each transmitted pulse 
actually resulted in three received pulses. The first pulse, 
which arrived almost immediately, was the direct source/ 
receiver path. The second was the reflection from the plate. 
while the third was the desired diffraction from the edge. 
Unfortunately, the first two pulses were considerably higher 
maeamplitude than the diffracted pulse. Even with the ability 
to exclude the first two pulses in the time domain, the 
ringing of the reflected pulse interfered with the analysis 


mmmeie Gittracted pulse. This limited the geometry to 
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Bemrce/receiver locations near the plate, which presented no 
Meer cular problem since this was the region of most interest 
anyway. 

Mmiteesnesults of this experiment are presented in Figure 16, 
math the free-field pressure at 2r used as the reference. 
@emcne COMparison, the theoretical B-T and H-K results are 
Shown for angles approximating those of the measurements. 
(These theoretical results are the same as those presented 
mieereures 5 and 6, now plotted as a function of frequency 
for specific angles.) The B-T theoretical results are 
plotted as a single line since they vary by less than one 
Memmemer the range of angles presented. It is clear that the 
measured data tend to agree with the results predicted by the 
Meeesolution rather than going to zero at the surface of the 
feemieomas in the H-K solution. The considerable scatter in 
mmemimeasured data, especially near 9 kHz at 30 degrees, is 
Mameroitcable to interference between the diffraction pulse 
@memeeiye ringing of the reflected pulse. 

ite measured data not only nee Ehcmp- 1 SOlUtION but 
meeomcontirm the presence of a scattered wave near the sur- 
mace. lf this scattered wave also exists on the surface, 
as predicted by B-T, (there is no reason to think otherwise) 
Meepoimts Out a serious fallacy in the Kirchhoff assumption. 
mmeorder to evaluate the integral in Equation II -25, the 
Wave field at all points on the surface of integration had 


to be known. This was accomplished, through the Kirchhoff 
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Mesumption, by simply taking the wave field on the surface 
membe the reflection of the incoming wave. This is equiva- 
lent to a boundary condition which forces the diffracted 
Wave field due to the edge to zero at the surface of the 
Meece When the solution is transformed to the time domain 
Meduations III-6 and III-7) where the diffracted and re- 
flected terms are considered separately, the diffracted 
term must go to zero at @ = Q° in order to satisfy the 
mm@esed boundary condition. That this is clearly not the 
meen indicates that the H-K formulation should not be used 
to estimate the backscattered wave field near the surface 
Maeeme Dlate. For forward scattering, the H-K diffraction 
Semue1on would also be expected to go to zero near the plate. 
Since Bremhorst's data [Ref. 12] show the presence of a 
forward scattered wave at the surface, the above argument 


Sameprobably be applied to this case as well. 


65 





feet NE TE PLATE, FORWARD SCATTERING EXPERIMENTAL RESULTS 


fee GENERAL 

When a rigorous diffraction solution such as Biot-Tolstoy 
tS applied to the problem of diffraction by a real object, 
one must be careful to ensure that the original boundary 
Bemere1Ons are applicable. For a finite barrier, two impor- 
tant considerations are whether the barrier is large enough 
jmembe cCOnSidered infinite in extent and thin enough to ignore 
thickness (the double edge). Since the geometry of the half- 
meme satttraction problem is such that everything is described 
an terms of a cylindrical coordinate system with the origin 
at the edge, there is a tendency to think of diffraction as 
mepeyean edge effect. This is the interpretation of the 
H-K boundary diffraction wave. This is also the simplist 
mmeerpretation of a time domain, impulse solution such as 
B-T, where it is easy to visualize an expanding spherical 
wavefront intersecting the edge and reradiating acoustic 
mmerey a5 it "propagates" along the edge. If this were 
generally true, one could solve any diffraction problem by 
Meseribing the diffracting object in terms of a series of 
edges (or, more appropriately, wedges) and adding their con- 
tributions. The problem with this approach is that the 
original solution assumed the presence of a rigid surface 


at some angles 8 = 0 and 8 = 8 extending to infinity in the 
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radial and axial directions and an infinite fluid elsewhere 
Beeeerigure 1). When the physical situation approximates this 
Bemdition, the diffraction problem can generally be treated 

by considering a source at the edge. For an extreme example 
where this is not true, one could consider the difference 
memween the diffraction due to a rigid half-plane and that 

due to a large plate shrunk to become a thin wire at the 
Tocation of the edge. 

A related problem concerns the treatment of barrier cor- 
meee tere the barrier may be "high'"' enough to consider the 
radial dimension to be infinite but its width (the axial, or 
@eeammension) is finite. A roadside billboard or a tall 
feeding Could be described this way. If the ground is 
totally absorbant (no reflections) Medwin [Ref. 8] treats 
miemaititraction from this type of a barrier by adding the 
impulse response of each of the three edges, truncating the 
solution for each edge when the wave reaches the corner, to 
obtain the total time domain solution. The assumption 
momcit here too is that the boundary conditions of the 
theoretical B-T solution are sufficiently well approximated 
By the actual barrier. 

The third major difference between the theoretical B-T 
feeution and the actual barrier is the double edge caused by 
the finite thickness. The theoretical solution is based on 
a rigid, infinitely thin screen whereas the actual barrier 


is often approximately rigid but of finite thickness. 
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meme researchers have minimized this problem in their mea- 
emements by placing a "knife edge" at the top of the barrier 
Weft. Kawai, Ref. 11] but this introduces additional sources 
Smee raction where the face of the barrier changes angle. 
In any case, the real barrier will always have some finite 
mirexness dimension and it is desirable to be able to calcu- 
late or at least understand the nature of its effect on the 
Mc iracted field. 

Pesceriles Of experiments was performed in order to define 


m@ese ‘finite barrier'' effects. 


Peer ORWARD DIFFRACTION AT A CORNER 

Bremhorst [Ref. 12] showed excellent agreement over a 
Wide frequency range between the B-T infinite plate solution 
Memmernise measurements on a 3/16" thick plate, provided the 
mepese) 1S Sampled in such a way as to exclude late arriving 
ditfracted pulses from other edges. However, discrepencies 
between measurements and theory of 5-10 dB were observed at 
frequencies approaching 50 kHz with source and receiver near 
the plate. The purpose of this experiment was to measure 
Mmmiemodittracted pressure field in the vicinity of the corner 
mae olate and compare it to the field predicted by B-T 
M@eane Medwin's truncation approach. The present geometry 
Was chosen to ensure that only the edges intersecting at the 
mormer of interest would contribute to the measurements 


(Figure 17). 
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The source was placed at the coordinates (X. > y : 


S 
fee. > cm, 8.0 cm, 8.0 cm, resnectively. Separate measure- 
memrs were made using the 4 and 1 inch sources to achieve 

the widest possible frequency coverage. The receiver was 
Meertioned along the x axis c= = 0,2, SsUjeat) Various points 
Memween r = 3 cm and r = 45 cm. This geometry was dictated 
by two considerations; source/receiver directivity and co- 
Mement Interference. Source and receiver are both to be 
Smart rectional in the theoretical B-T solution. In order 

tO approximate this adequately, it was considered important 
that the source appear omnidirectional when viewed from 
points along a significant portion of the edge near the 
corner. By orienting the source and receiver in this way 

all points on the edge between the corner and the point where 
the least time path crosses the edge were within 2 dB of 
omnidirectionality up to approximately 30 kHz. The symmetry 
of the geometry was chosen to eliminate the frequency domain 
interference pattern created by taking the Fourier transform 
of a time signal which has two coherent pulses separated by 

a time delay. The geometry ensures that the diffracted pulses 
from both edges arrive at the same time. This was accom- 
plished by moving the receiver slightly to maximize the total 
pulse signal at the nearest and farthest ranges and locating 
the intermediate measurement points on a straight line 


between them. 
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These considerations resulted in a geometry which was 
E@mewhat artificial but would allcw a relatively simple analy- 
Siemeot the corner effect. The Fourier transformed diffrac- 
mea pulse data for the 1"' and '"" sources are presented in 
Figures 18 and 19 respectively. The free-field pressure at 
menranee Of 28 cm (the distance from source to corner) was 
used as a reference pressure. Although the measured data 
using both sources show good agreement in the frequency 
region where they overlap, their dimensions are large enough 
to be considered as an effect on the results and are there- 
mmmempresented separately. The theoretical results from 
truncating and adding the B-T solutions for the two 
Mem@oecndicular edges is also shown. The agreement between the 
measurements and this superposition of the truncated 
meumm@oments Of two cases of B-T theory is seen to be excellent 
jemoss the entire frequency region. 

The theoretical and measured results in Figures 18 and 
meealso Clearly show the effect of moving closer to the 
shadow boundary. Because of the geometry shown in Figure 17, 
decreasing r also moves the receiver nearer to the shadow 
boundary. The theoretical solution as the receiver approaches 
mime shadow, is shown to approach the delta function form 
(Appendix A). The measured data clearly show this result, 
meethey approach a flat spectrum for small r, with the more 
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mamiliar f frequency dependence observed at larger r. 
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In Figure 20, these same measurements are presented as 
[ereracted pressure as a function of range (from the corner) 
Meupeese lected frequencies, once again compared to the sum of 
menumeated' B-T theoretical solutions. Also presented for 
comparison are lines of arbitrary level with slopes of -10 
meer. cylindrical) and -20 log r (spherical). An apparent 
molmece eat the corner would be expected to radiate spherically 
fee an infinitely long edge would radiate cylindrically. 
Because of the geometry the ranges to the corner and to the 
meme On the edge where the least time "ray" crosses are 
essentially equal so the data may be examined directly for 
mem@ience Of a spherical divergence from the corner. 

Mmemaata, as well as theoretical results, of Figure 20 
range from a slope of less than 3 dB per double distance 
Mear the corner, low frequency, to almost 6 dB per double 
Meeeeance at the highest frequency and largest distance from 
mmemmecorner. in terms of the dimensionless parameter kr 
Miavenumber times range), it appears that the corner cannot 
Depcmaracterized as a simple spherical or cylindrical source 
fmeerow Values of kr. As kr is increase, the range depen- 
Menee gradually increases to avalue approaching 6 dB per 
double distance. Although the results are somewhat incon- 
mesive due to limited data, the two line sources, one at 
Gach edge, begin to resemble a virtual point source at 


higher values of kr. 
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An additional experiment was performed to better define 
mememrterence between diffraction by an edge and diffraction 
Dy a corner. The plate used in this experiment was identical 
memeemae used for the corner (Fig. 17) except that the source 
Gmemreceivyer were set up on a long edge, as shown in Figure 
Biers before, the measured data are compared to the B-T 
@ebeulated results (for a single infinite edge) in Figures 
22 and 23. Unfortunately, the agreement between measurements 
m@emuieory iS not nearly as good in this case as it was for 
the corner. The reason for the significant deviation between 
measurement and theory at higher frequencies and small ranges 
1s not known. Figure 24 presents a direct comparison between 
memmer and infinite edge results, as a function of range, 
along with the appropriate B-T theoretical results. 

From this comparison it can be seen that the diffracted 
pressure field behind the corner is approximately 6 dB higher 
mieamecnat behind a half-plane (single edge) for the larger 
ranges and higher frequencies. As the receiver moves closer 
to the plate, the difference between corner and single-edge 
data decreases. Furthermore, there does appear to be a fre- 
quency dependence in the difference between the corner and 
edge at a given range. To better define this frequency 
dependence and to estimate the effect of moving away from 
the shadow boundary, a final experiment was performed. 

Referring to the geometry of Figure 17, the receiver was 


Miaced distances ({r) of 18 and 32 cm from the corner on a 
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imememperpendicular to the plate (this is the same geometry as 
the previous corner measurements), Measurements were also made 
mememoeances Of 18 and 32 cm from the corner on a 45° line 
Sepemaing into the shadow region "behind" the plate as shown 
in Figure 25. As before, four identical measurements were 
feeematoneg a Single "infinite'' edge for comparison. The 
(eomeeuae of the Fourier transformed data from the corner, 
Givided by the corresponding data from the single edge, are 
Meeemmeed in Figure 26. The upper two curves present data 
woeemmon the perpendicular line while the lower two are taken 
on the 45° line (deeper in the shadow). Also shown with the 
[ieemeeurves iS the theoretical difference between the corner 
omemeee@oce based on B-T. The general trend indicated by the 
MmeemeewoO CUrves is one of less than a doubling of pressure 
(due to line sources) at the lower frequencies, increasing 
Eememementiy more than double the pressure at the higher 
frequencies. 

When the receiver is moved deeper into the shadow region, 
Mmiemenaracter of the ratio plotted in Figure 26 changes 
dramatically, as shown by the lower two curves. The monotonic 
frequency dependence seen near the shadow boundary has been 
replaced by a more complicated dependence. The average value 
eeeeme difference appears to be approximately 6 dB, even at 
mem@enlower frequencies. The detailed frequency dependence 
will be addressed in the next section. The general conclu- 


pon to be drawn from Figure 26 is that the farther the 


81 





Se 
Near Shadow, 
Boundary = 
Receiver 


Deep Shadow 
(Gragonal } 
Receiver 





° ’ srl! 
mreure 25. Experiment Geometry for "Near Shadow Boundary! 
"Deep Shadow" Forward Diffraction Comparison 
Measurements. 


82 





ATT 


NEAR SHADOW BOUNDARY (Fige 25) 








== Vel erm 
= 
‘4 
a 
oS 
oy 
Ou 
r>3a em 
+ ¥ 
as Truncated B-T Theory 
DEEP SHADOW DIAGONAL (Fig. 25) 
+ 3 + 
ot: 
ee tt Le 
T7219 cm 
- 
OS 
S 
- 
3 
. c 23a cm 





a = 10. OP 
Frequency in kHz 


Figure 26. Comparison of P__, as a Function of Frequency 
for Receiver LocStions in the Deeo Shadow and 
Near the Shadow Boundary (see text). 





memeeiver moves into the shadow region, the more the diffrac- 
ted pressure field looks like the coherent addition of two 
infinite edges, plus a small freauency dependent variation. 

In summary, the measurements presented herein show that 
Ewemeasic approach of adding impulse diffractions from each 
edge in the time domain is valid. Near the corner, the 
diffracted vressure is less than a simple addition of two 
infinite edge diffractions. In this region, the apnroach 
proposed by Medwin [Ref. 8], wherein the time domain solu- 
tion is ''truncated" at the point where the edge wave reaches 
the corner, estimates the diffracted pressure verv closely 
as shown by Figures 18 and 19. The usual treatment of the 
COrmer as a source of spherically diverging diffracted sound 
1s not evident in the measured or B-T theoretical results, 
at least over the lower range of kr presented. At high 
values of kr, the results mav be approaching spherical 
Pee roence. 

In the high frequency limit, Keller [Ref. 16] predicts 
that the corner does produce spherically diverging waves but 
that the resulting diffracted field decreases with increasing 
EPrequency much more rapidly than the field produced by an 
edge. This characteristic has been used as justification 
memeaisregarding the diffraction at a corner when approaching 
Practical problems [cf. Ref. 17]. Figure 26 would seem to 
Indicate that radial distance from the corner and frequency 


meeemot Sufficient to determine when to disregard the effects 
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Sema cOrner; one must also consider the relationship of the 
Meceiver location to the shadow boundary. The dimensionless 
M@apameter implied in Medwin's truncation anproach is kd where 
@emsecne distance from the corner to the point where the 
Beast time "ray" intersects the edge. No substantiation 

for this parameter is offered here, however it does take 
into account both distance from the corner and nearness to 


the shadow boundary. 


See ORWARD DIFFRACTION BY A STRIP 

The next experiment was designed to investigate the be- 
menor Of diffraction from an ‘'infinite" strip of material 
whose width was of the order of a wavelength. Two strips 
Were actually measured, with widths 10 cm and 4 cm. Given 
S@emmeasurement capability of 1.5 to 40 kHz, these two strips 
provided a total range of ka (a represents the width of the 
memo Or 1.09 to 72.8. Source and receiver were located on 
Segerpendicular line through the center of each strip, as 
shown in Figure 27a. The source was 14.5 cm from the strip 
for all tests while the receiver was located at distances 
meee Strip of 9.5, 235.2, 39.5, and 59.5 cm in the case 
mie 10 cm wide strip and 9.5, 23.2 for the 4 cm strip. 
Once again, this geometry was chosen to eliminate the cohe- 
mene interference that would result from different path 
Meneths around each edge. Identical measurements were also 


Mertormed on a single “infinite'' edge as shown in Figure 2/b. 
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meeure 27. Strin Diffraction Experimental Geometry. 





As long as linear superposition of edge effects holds, the 
diffracted pressure at any of the measurement points will be 
double for the strip compared to the corresponding plate 
(single edge), since there is no truncation involved. With 


Situeein mind, all data are presented as the following ratio: 


P Brcimuenacved pressune behind strip 
ea t diffracted pressure behind plate 


When superposition holds, there will be a simple doubling of 
Pressure and this ratio will be +6dB at all freauencies. 

This ratio, calculated from the measurements on the 10 
cm strip and corresponding plate, is presented in Figure 28. 
iitemcdata are not presented above 30 kHz because of limita- 
tions in the determination of the perpendicular to the strip 
(i.e. the path lengths around the two edges could not be 
der zed well enough to vrevent cancellation effects above 
this frequency). The oscillating pattern in the low frequen- 
meportion of the measurements is evident from Figure 28. 
memes 1S interpreted as coherent interference between two 
Petals, two features can be extracted eer the data. First, 
this interference pattern does not depend on distance from 
mmeestrip/plate. Second, by estimating the spacing (in fre- 
quency) between adjacent peaks or between adjacent nulls the 
Meearent path difference can be calculated from the relation 
d=} = C/ af ime ceus wf bom eleCe@dedmor rigure 238. d.15 


estimated to be 9 cm, quite close to the strip width of 10 
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M@eeure ¢O. Comparison of P as a Function of Freauency 


at Various Disthfkes (ro) men the 10 cm Strin 
(Geometry of Fig. 27). 
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Mammo Verity that strip width is the characteristic dimen- 
emt, the experiment was repeated for the 4 cm strip. 

fies + CM strip results are oresented in Figure 29. The 
Semepewidth from Figure 29 data is estimated to be 4.1 cm, 
Semcdiming that the width of the strip is the factor deter- 
ieniine the form of the “interference pattern". 

mieedata measured at 9.5 and 23.2 cm are averaged and 
Meeeead 2S a function of ka for both strips in Figure 30. 
Meeeeemiplicity, the data from both 1 inch and 1/2 inch sources 
were averaged in the overlapping region 6-10 kHz. Figure 30 
clearly shows that the frequency dependence is directly re- 
[Mieed to the strip width a, with a ka spacing of 2n between 
peaks as expected from a simple interference pattern. 

The phvsical explanation offered for this interference 
aS Called secondary diffraction or secondary scattering. 
mmememp lication of secondary diffraction to the strip is 
memuistrated in Figure 31. The figure shows only the 
interaction at one edge; the mirror image of the process 
shown 1S occurring simultaneously at the other edge. The 
feeeramlines in Figure 31 depict the path of the pulse as it 
diffracts over the upper edge, while the dashed lines show 
how the diffraction of the same pulse at the opnosite edge 
generates signals on both sides of the strip which propagate 
Memeeoethe Strip and diffract once again at the upper edge. 
imeedoubly diffracted pulse arrives at the receiver at time 


equal to the width (a) divided by the propagation velocity 
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Figure 29. Comparison of P 


Some Me cuOn Of ss requency 
at Various Distahées Go cor the ecm Strip 
(Geometry of Fig. 27). 
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Meyeeeacer than the initial pulse, resulting in the observed 
[Mieemeerence pattern. This late pulse would have two dis- 
mu@t@eemcenaracteristics, arising from the additional diffrac- 
tion, which are immediately observable in the measured data 
Seeeure 30. First, it will be considerably reduced in 
ampemtude relative to the earlier (single diffraction) 
pulse. This is seen in the data as a weak interference 
pattern (the peaks are approximately 6 + 1 dB compared to 

6 + 6 dB if the two were of equal amplitude). Second, the 
difference in amplitude between singly and doubly diffracted 
pulses will be greater with increasing frequency since each 
diffraction should increase the freauency dependence and 
decrease the strength of the diffraction by approximately 
fo This causes the decaying amplitude of the interference 
pattern with increasing frequency. 

A method of calculating the secondary diffraction has 
been provosed by Medwin? and imelemented by Ms. Emily Chiids. 
The general description of the implementation which follows 
1s due to rete When a spherical pulse intersects an 
edge, a boundary or edge wave is generated which anpears 
to propagate along the edge. Calculating the B-T time domain 
solutionin discrete time increments, this edge wave becomes 


a series of discrete sources on the edge, each with a 


Meer sonal communication to author from Dr. H. Medwin. 
Personal communication to author from Ms. E. Childs. 





Sa@eacteristic amplitude and phase. Extending this to a two 
BPdce System (the thin edged strip in this case), these dis- 
eememsources on the first edge (now called secondary sources) 
Stemearitract at the second edge and contribute to the dif- 
fracted field at the receiver. The proner amplitude is 
meemeomed tO each secondary source through a system of vir- 
tual receivers located at the second edge. The phase of 
Semmeeource iS taken into account by performing the entire 
calculation in the time domain. Figure 32 presents the 
Geers Of this theoretical calculation in terms of Deaee eng 
Stemestrip at a receiver distance of 25.2 cm., compared to 
the measured data from Figure 30. Two theoretical curves 
are presented for each striv; one which considers the 
woeomedry diffraction on both sides of the strip (double 
secondary path) and one which considers secondary diffraction 
only on one side (single secondarv path). Although the 
double secondary path is the intuitively correct 
interpretation, the measured data obviously agree with the 
Single path approach. One possible explanation is that in 
uSing the double path interpretation, each path should be 
reduced in amplitude bv 1/2 (the diffracted pressures along 
the paths on either side of the strip are of equal amplitude 
but opposite sign). Keller [{®ef. 16] describes a similar 
phenomenon for the thin slit in an opaque plane using the 


meemeerical Theory of Diffraction. 
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fyeeenis point it is also appropriate to point out that 
emestse of secondary diffraction in the preceding explanation 
Beene interference pattern would not be possible with the 
mime ltz-Kirchhoff approach, since the scattered or dif- 
fracted field on the plate would be zero. There would be no 
way the pulse could propagate between edges. 

iemappears that the secondary diffraction explanation can 
also be applied to the somewhat more complicated geometry of 
@@emeorner. Recalling Figure 26, there was a detailed fre- 
quency dependence in the results, especially where the 
receiver was deep in the shadow region (the lower curves of 
Figure 26). Consider now a corner with one edge horizontal 
and one vertical as shown in Figure 17. There will be a pulse 
tEMat travels from source to receiver by diffracting over the 
iemizemctal edge. In addition, a pulse also diffracts at the 
Memtrcal edge, propagates across the plate near the corner, 
am@a@madittfracts over the horizontal edge to the receiver. 
Mmeuke the strip, the difference in path length of the two 
mmeseoeis a function of distance from the corner, an effect 
moecnmecan be observed as a difference in the interference 


Meeecrns in the two lower curves of Figure 26. 


D. FORWARD DIFFRACTIONS BY A THICK BARRIER 
iiesusual theoretical approach to solving the problem 


Memeeterraction by a plate is to consider it to be an 
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Meemitely thin, rigid half-plane. However, a real barrier 
Mmmeaesitenificantly attenuates the direct path of sound must, 
Geeeecessity, have a thickness which is a significant frac- 
tion of a wavelength at high frequency. In Reference 12, 
Bremhorst presented data which showed the diffracted pressure 
at various points in the shadow of a 53/16" thick steel plate. 
In particular, Figures 37 through 46 of Reference 12 showed 
that the measured diffracted pressure was generally less 
than that predicted by B-T, especially near the surface of 
the plate (on the shadow side). Bremhorst attributed this 
Mommie tinite thickness of the plate. To confirm this 
hypothesis, an experiment was performed using two barriers 
of different thickness. 

The variable thickness was achieved by placing a piece 
of aluminum stock of the appropriate thickness between two 
3/16"' aluminum plates and machining the resulting "edge" 
smooth and square. Damping material was placed in the gan 
between the plates to attenuate the direct path. An example 
of the thick edge can be seen in the photograph of Figure 
14. Thicknesses of 2.2 and 3.5 cm were evaluated utilizing 
the geometry shown in Figure 35. Measurements were made at 
8 = 300° and 345°, where there was considerable disagreement 
between Bremhorst's measurements and B-T theory. The measured 
meereare shown in Figure 34, compared to the Bremhorst's 


mesults and the B-T theorv. 
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Figure 33. Thick Barrier Measurement Geometry. 
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Pressure in dB re free=field at 50 cm 


Measured Data: 


3d -£5 48 cm (Bremhorst) 
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Figure 354. Forward Diffracted Pressure as a Function of 


Frequency for Various Thickness (t) Plates, 
aS ce Pen o en. 





Miesdata presented in Figure 34 indicate that barrier 
thickness 1S an important parameter in determining the dif- 
meaeted Nressure field, especially when the thickness is on 
fiemorder of one wavelength (10-15 kHz for the thick barriers 
tested). It has been suggested by Medwin : fale Set Ad CK = 
Mesoecan be handled by a secondary diffraction calculation 
meiseason the B-T theory, as discussed earlier for the thin 
strip, by considering the thickness to be a double edge. 
mesmias not been verified at this time, although the thick- 


ness does appear to explain the discrepancy between Bremhorst's 


data and the B-T theory. 


E. FORWARD DIFFRACTION BY A BARRIER ON A RIGID BASE 

Peetyoical noise control barrier consists of a plane 
screen or plate mounted perpendicular to a flat base, as 
smermeim Figure 355. The total pressure field at the receiver, 
Hiwit 1S in the shadow region, consists of the source signal 
arriving via multiple paths. The usual approach to solving 
Ears problem is to consider that for the top edge there are 
four geometrically distinct paths, each involving one 
Serraction at the top of the barrier, made up of all vossible 


combinations of source, receiver, and images of source and 


Personal communication to the author by Dr. H. Medwin. 
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Meeemver in the base, as shown in Figure 35. Using the no- 
@eenon Of Figure 35, these paths are denoted as SDR, S'DR, 
freee DR’. It should be pointed out that S' and R®' here 
ieemenue images as used in optics and acoustics to describe 
jmeebection from a plane surface. The term "image" is also 
meeamoy authors in the noise control field [cf. Ref. 5] in 
Eonjunction with the "image of the source in the barrier", 
as discussed in section II. This usage is considered here 
to be inappropriate and will always be set apart bv quotation 
marks ("image'') to prevent confusion with the more common 
term (image). 

iomcest the hypothesis that the four paths described 
above adequately characterize propagation from source to 
mememver in this case, a simple experiment was performed. 
Pweparrier with base was constructed of 3/16 inch aluminum 
plate with an overall width of 120cm and barrier height of 
mepemean the basic configuration shown in Figure 35. Source 
and receiver (both 1/2'' B+K microphones) were positioned so 
meeemecacn of the four paths would have a different length, 
Miameach path would intersect the axis of the microphone 
@eeapproximately the same angle. This would allow at least 
megialitative comparison of the relative strengths of the 
mirdews paths. Use of the pulsed source of earlier experi- 
ments would allow the identification, in the time domain, 


meweach path by its characteristic propagation time. 





Samiiarly, undesirable paths such as backscatter from the 
mieeeor the finite base and diffraction around barrier sides 
would have pronagation times longer than the longest (S'DR') 
Sarieor interest and could be eliminated in the time domain. 
The results of this experiment are shown in Figure 36, 
meomgteed tO the calculated arrival times of each of the four 
Mieioe Lt is clear from these results that the four paths 
constructed from the source, receiver and their images in 
miemeace are the dominant vaths in this basic barrier confi- 
eeteron. iif any other paths exist, they must be of rela- 
tively low level or coincide with one or more of the four 
Paomeswpaths. The path due to the "image of the source in 
the barrier'' as discussed by Isei [Ref. 5} would coincide 
With the SDR and SDR' paths of Figure 36. The relative 
amplitude of the various paths also behave as expected, with 
the later arrivals (longer paths) having lower amplitudes. 
The strengths of the various vaths were examined bv 
varying the geometry slightly. Bv increasing the barrier 
moemeeies to 42 cm and changing the receiver location (Figure 
faleeche Separation in time between the paths was increased 
Emeouen tO allow a separate spectrum analysis of each pulse. 
The data of Figure 58, taken with the source/receiver oriented 
to insonify/receive both SDR and S'DR paths equally 
(approximately), show the separation between the two pulses, 
macdthe Seneral trend in their amplitudes. As expected, the 


S'DR pulse appears to be an exact replica of the SDR pulse. 
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Mmemaretiual strength of each path was measured hy "pointing" 
Meemoource in the direction of the path of interest (the re- 
Seaver always pointed toward the edge, D) and Fourier trans- 
moummne Only the received pulse corresponding to that path. 
itmoewas done for both the SDR and S'DR paths. The data, 
normalized by the free field pressure at the receiver with 
memeerrier and base removed, (source/receiver pointed at 
Sameormncr) are presented in Figure 39. Also presented for 
weeeeisoOn are the B-T infinite half-plane theoretical results 
for each path. The agreement between theory and measurements 
Soermeeim Figure 39 is excellent, indicating that the diffrac- 
meemerreld at the receiver due to SDR and S')DR paths may be 
Gomouced separately, using the geometry of the source and the 
image of the source in the base. The total diffracted 
pressure at the receiver may then be obtained by adding the 
Semertoucions of each path. 

ments important to note that the theoretical results in 
Figure 39 were obtained from the complete form of the B-T 
meedte half-plane solution. The complete theoretical 
mmmmenron tor the diffracted field at R in Figure 35 consists 
Bemtne sum of four terms. Isei [Ref. 5] approaches the same 
maepbem by adding the contributions of six different paths, 
Maesitown in Figure 3 of Reference 5. It anpears that he 
has used the solution due to Macdonald in a form which is 
Beparated into "real" and "image" components as discussed in 


section II and is applying these separately to the six paths. 
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oe Ove AO 5G. 
Frequency in kHz 
Diffracted Pressure as a Function of Frequency 


for the SDR and S'DR Paths Compared to the 
B-f Theory (Geometry of Fig. 37). 
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ieee 59 Shows that the correct approach is to applv the 
Semobete SOlution (e.g. in Macdonald's solution this would 

be the sum of the "real" and "image" source components) to 
Seamer the four paths to obtain the correct diffracted field. 
Siemmeed not consider the "image of the source in the barrier" 
Mem@awse it 1S an integral part of the complete solution. 

Meaddition to the four paths described earlier, it is 
possible to have an additional path which propagates from 
Mmeemesource, diffracts at the intersection of the barrier and 
fieere>, then diffracts over the top edge of the barrier to 
memereceiver, shown in Figure 40 as SBDR. This assumes a 
Mompeeerpendicular intersection of barrier and base since 
emeeeoy {Ref. 7] shows that there will be no diffraction from 
a 90° "interior" corner. The theoretical prediction of this 
Meeeeracted pressure at the receiver was approached in the 
peemway as for the thin strip, discussed earlier. 

The geometry shown in Figure 40 was constructed of 3/16" 
aluminum, as before, with the source and receiver located as 
Peeve fhe SDR and SBDR paths (see Figure 40) can be seen 
directly in the time domain in Figure 41. (The relative 
amplitudes are only approximate due to source directivity.) 
Once again, the source was pointed toward the path of interest 
and the received pulse was analyzed separately. The frequency 
meeetrum of the SBDR nulse, normalized by the direct, free- 


meld pressure, is shown in Figure 42. Also shown is the 
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Figure 41. Received Time Domain Signal Due to the SDR 
and SBDR Paths (Geometry of Fig. 40). 
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mieeretical result from the secondary diffraction program. 
femecomparison, the SDR path data and theorv from Figure 359 
meemarscO presented. (The SDR geometry is identical to the 
memes experiment.) The additional diffraction in the 
eeweaepath is seen in the data as a different frequency depen- 
@ee@ec as well as a lower overall level. When the B-T theory 
memocad 2S a basis for secondary diffraction, there is 


Beweetlent agreement with the measured data. 


dike 
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Meemexact impulse diffraction solution of Biot and 
iescoy (B-T) has been compared to the exact cy  diffrac- 
iemeoOolution of Macdonald for forward scattering and to 
the impulse solution to the Helmholtz equation (H-K) for 
Pemeseatter using the Kirchhoff approximation. For back- 
Seeeeeer from an infinite edge, B-T and H-K have been found 
momemeeance to a Similar although analytically different time 
fee form. Comparison of B-T and H-K in the frequency 
domain showed substantial disagreement, except for geometries 
feeteemeune diffracted pulse was close to the reflected pulse 
in time, where both solutions approached the same value. 
Peum@ieseatter near the plate was found to be grossly under- 
estimated by the H-K solution, due to the boundary condition 
Beeamoy Kirchhnoffto simplify the integral. The backscattered 
pressure near the surface predicted by B-T was confirmed by 
laboratory measurement. 

macaonald and B-T solutions for forward diffraction from 
a plate were both found to be in good agreement with the 
meeemament. Both solutions can be expressed as the sum of 
a term containing the source/receiver angles as (9-8) and a 
meemecontaining (9+). These terms were originally identi- 
fied by Macdonald as the "real't and "image" source terms. 


Mmemeoncept of "image of the source in the barrier", as used 


114 





Meme NOise control literature, is found to be a mis- 
Mmieenpretation of Macdonald's solution. For a barrier on 
eeeaerd base, the diffracted pressure due to the real source 
and the true image source (the image of the source in the 
base) were measured separately and both were found to agree 
with the complete (both ae ad eo temic) esa SO LutTom. 
mire solution of finite barrier problems was found to 
be facilitated by the time domain, impulse nature of the B-T 
fommeron. Measurements of forward diffraction near a plate 
corner show that Medwin's proposed approach of truncating 
Siemo-I time domain solution to account for the corner accu- 
Mmemety predicts the diffracted field. Experimental results 
Bemm@e agree with simple spherical radiation from the corner 
Srememost Of the range of kr between 1 and 150. Limited 
data do show a trend toward spherical divergence as kr 
@eeeoadenes 150. Measurements of the forward diffraction by 
fem Strip confirm the presence of secondary diffraction, 
imeemmewaves that first diffract at one edge, propagate across 
miemsurtace of the strip, and diffract again at the other 
edge. The time domain B-T solution is shown to be a useful 
Meweeing point in characterizing this secondary diffraction. 
At high frequencies, forward diffraction loss by a plate of 
memte thickness is found to be considerably greater than 
mieetheoretical prediction, for an infinitely thin screen. 
It appears that this may also be explained by the double 
edge creating secondary diffraction with an additional 


Peeetraction contribution at the second edge. 
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merROXIMATION OF THE BIOT-TOLSTOY SOLUTION 
NEAR THE SHADOW BOUNDARY 


peepee the B-T solution is a delta function impulse solu- 
tion, it 1s expected that it will approach the delta function 
form for limiting cases where the receiver is approaching 
the boundary between "illuminated" and ‘shadow" regions, 
herein referred to as the shadow boundarv. As this limit is 
@eproached, the early time portion of the diffracted pulse 
will increasingly dominate the solution, creating a difficult 
Mumerical calculation problem. To better understand this 


problem, the solution will be approximated for the following 


Pomentions: 
A. ar oem thin plate) 
reeeeo = 6 + T + € (receiver in the Shadow region, ¢e 
: small) 
fee - t +0 from the positive side (early 
: diffracted pulse) 
BOT a = 27, the complete solution can be written as 
ot 
p(t) = - 52° {e}[{rr_ sinhY] [exp(-'/2)] (A-1) 
Sans _ 
where 
Dos 2 2 ee 
io. = arg Gosh poet = eee oe (A-2) 


27k 
O 


and 8 is a four term sum written symbolically as 
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i 
sin tp (moto) | 


(BRB 
i} 


1 - 2exp(-5Y)cos[$(n#9#8,)] + exp(-Y) Eee) 
mene the early diffracted pulse 

oee- t - tT +0 

0 


Medwin [Ref. 8] shows that when Eccl, Cle sEderorucOntain ing 


0 
Sinh Y can be approximated as 


tr 


-1 

oT Eoin) = 2 clea 
Meeisesame result for the specific case of backscatter was 
eerie in section III.} This factor goes to infinity as At 
approaches zero and is generally assumed to be the factor 
that determines the time dependence of p(t) under these 
Circumstances. Time dependence also enters A-1l through the 
exp(-Y/2) term but this will be ignored because this function 
only varies between 1 and zero and is relatively constant 
for Y (and therefore At) approaching zero. 

imme 8 term contains both time and geometry dependence in 
feeeener complex wavy. The approach here will be to separate 
the 8 into four terms and examine each one in the limit as 
both Y and € approach zero. (Y is directly, although not 
memedriy, related to At. The expansion of arg cosh is not 
trivial and so the limiting procedure here will be with 
meewect to Y rather than t.) The 8 term can be written as 
(8 er 8 + 8 oes) Withee nemeWlinamneEsmieterring sto the 
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Various combinations as shown in Table (A-1). 
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iene the usual, first order approximations for sine and 
cosine of small angles, these terms can be approximated by 
the following 


Sim 0. “ee / 2° cos. oO 
O 0 


ae. = ee fo. 
it) 2 exp /2)[cos8 . a ean oF.) 1 eco =) 

; Zz “f2 

i 1 + 2exp(-¥/2) + exp(-Y) (A-4) 

fo 

i 
ee 2 CX (+1 / 2) eee xD (-Y) 

, sin 9, + &/2 cos % 


1 - 2exp(-Y/2) (cos 8 oe eran 8) tone xp) 


With these first order approximations, the behavior of 


each 8 term can be examined at Y= J 


Meese oC cooe es 


on. | z 

a Y=0 Z2{1 + cos 6 - e/2 sin @ ) 

0 0 
g | » oo Se a geomet 
me Y=0 CES CO= Ce ete 
an | z- 5 
*-"y=0 : 
BS | = (denominator goes to zero) 
Y=0 


Memes first order analysis; ie iS seen that three 
of the four 8 terms have the potential of dominating the 
Meemesolution at small values of Y. These terms will now be 
examined to second order in the denominator. (The numerator 


contains no difference terms and so the first order 


ens 





approximation for small e« will be retained for simplicity.) 


Considering 8 _, the following approximations are used 


cos ie 72) 2 ee 








on 
Y2 
mmpe(-Y) = 1l1-Y +35 
momobtain 
o Ze —- 
2 Nea Woe 


— Ye 
i oS = 9 ao. 
Seweeeting terms higher than second order in the small 


Seamer ties Y and ec, 


Ze 


a ; 
eee (lot 2) — 


which approaches a finite value as Y approaches 0, 


im 8 , = - a 
Y+0 : 


Weane an identical approach yields the following results 


Poe, and 8 
; ey 
oa 2 Q 
lim 8,. = - am Save (A-6) 
oO ag . 
2(l+cos 6 + ~/8 cos 6. - ©/2sin @_) 
0 O 0 
Son 0 uae 2.COs 8 
im Ss = .— <=. ae 
0 2(1l-cos 8 es / Secs ae e/zZsin 8) (A-7) 
A detailed analysis of the 8,, and 8 __ terms was not 


done, however, they assume the following maximum values at 


m= (): 
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since these limiting values are the same as for the eee eran 
memwil! be less for other values of 8 > Meise SeOne liided. cia t 
the 8 , term will dominate 8 for the early diffracted pulse. 

For comparison, the (rr sinh one term can alse be 
fepeoximated for Y approaching zero as 


ae Y3 
(rr sinh ie = ey, - + ) og L 
O ae a 





(A-8) 


Comparing Equations A-8 and A-5 leads to the conclusion that 
the rr sinh Y factor will determine the time dependence of 
the early diffracted pulse only when e* is much greater than 
Semana Y iS much less than l. 

imeecerms of the basic diffraction problem geometry, this 
means that as the receiver approaches the shadow boundary 
from inside the shadow (e approaches zero), the time region 
over which the (rr sinh oe factor will determine the time 
Seemaence Of p(t) becomes very small. This region must be 
determined, based on the geometry, before an analytical plus 
@rettal approach such as outlined by Medwin [Ref. 8] can be 
employed. This 1s fliustrated in Preimecm.-15 4-2 when show the 
behavior of (sinh on AV Come ee oe SUG eG C Hane m aon a CO 
Memes Of €. Since the B-T solution contains the product of 
mmese terms, one can be assumed to be the factor which deter- 


mines the time dependence only when the slope of the other 
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Figure .—l Behavior of the E=(eSseolutuon tear the Shodow 
Orinoco. 6 = lows 
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Figure Ae2 Behavior of the.BeT Solution Near the Shadow 
boundaery,6- 4.2 . 
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mmemeaches zero. Based on this, the (sinh Y)"— term would 
dominate below approximately Y= .1 for e« = 18° and below 


Y=,01 fore = 4.2°. 
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